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Abstract 


This  report  contains  curves  and  formulas  for  the  calculation  of  the  critical 
external  pressure  of  finite  length,  circular  cylindrical  shells  with  sandwich 
walls.  The  facings  of  the  sandwich  are  equal  and  isotropic  and  their  indi¬ 
vidual  stiffness  is  not  taken  into  account.  The  sandwich  core  is  isotropic 
or  orthotropic  having  natural  axes  in  the  axial,  tangential,  and  radial  direc¬ 
tions  of  the  cylinder.  Curves  are  given  for  cylinders  with  sandwich  walls 
having  isotropic  cores  and  orthotropic  cores  having  certain  relative  elastic 
properties.  If  the  cores  are  very  rigid,  the  method  yields  results  that  are 
substantially  those  of  von  Mises. 


—This  Note  is  a  revision  of  a  report  by  the  same  title,  issued  May  1959  as 
Forest  Products  Laboratory  Report  No.  1869.  This  revision  was  made 
by  E.  W.  Kuenzi  and  B.  Bohannan. 
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—This  progress  report  is  one  of  a  series  (ANC-23,  Item  57-3)  prepared  and 
distributed  by  the  Forest  Products  Laboratory  under  U.  S.  Navy,  Bureau 
of  Naval  Weapons  Order  No.  19-61-8041  WEPS  and  U.  S.  Air  Force 
Contract  No.  DO  33(616)61-06.  Results  here  reported  are  preliminary 
and  may  be  revised  as  additional  data  become  available. 
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“Maintained  at  Madison,  Wis.  in  cooperation  with  the  University  of 
Wisconsin. 
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Introduction 


This  report  presents  design  curves  for  the  critical  external  radial  pressure 
of  circular  cylindrical  shells  with  sandwich  walls,  calculated  according  to 
the  formulas  developed  in  Forest  Products  Laboratory  Report  1844-B  (9)— 
The  sandwich  cylinder  walls  have  isotropic  facings  and  orthotropic  or  iso¬ 
tropic  cores.  The  natural  axes  of  the  orthotropic  cores  are  axial,  tangen¬ 
tial,  and  radial.  These  formulas  reduce  substantially  to  those  developed  by 
von  Mises  (7,  10,  14)  when  the  core  is  very  rigid. 


A  great  deal  of  investigative  work  has  been  done  on  isotropic  cylindrical 
shells  subjected  to  external  pressure  since  report  1844-B  giving  theoretical 
analysis  of  sandwich  cylinders  was  written.  It  has  been  found  that  experi¬ 
ment  sometimes  yields  critical  loads  that  are  less  than  those  predicted  by 
von  Mises'  theory  (13).  This  has  been  attributed  to  two  causes.  First,  the 
experimental  cylinders  contained  imperfections  that  lowered  the  critical 
load  (1,2,  3,8^,  12).  Second,  lower  energy  levels  are  associated  with  post- 
buckling  configurations  of  the  cylinder  than  with  those  just  at  buckling,  the 
former  being  reached  without  the  necessity  of  passing  through  (snap -through) 
the  latter  or  the  energy  necessary  for  passing  through  the  latter  being  supplied 
by  vibration  or  shocks  (4,  5,6,8). 


The  curves  published  in  this  report  do  not  consider  snap-through  buckling  or 
cylinders  with  imperfections.  Sandwich  cylinders,  however,  are  much  more 
perfect  than  their  solid  counterparts  because  they  are  thicker  and  the  effect 
of  an  imperfection  is  in  proportion  to  the  ratio  of  its  amplitude  to  the  thick¬ 
ness  of  the  cylindrical  shell.  Also,  the  curves  neglect  the  stiffnesses  of  the 
individual  facings.  These  stiffnesses  add  to  the  critical  loads  when  the 
cylinders  are  short,  and  it  is  for  short  cylinders  that  "snap-through"  is 
likely  to  occur  (6). 


Development  of  Formula  from  Which 
Design  Curves  Were  Calculated 


The  critical  external  pressure  is  found  by  solving  equation  (51)  on  page  23  of 
Report  No.  1844-B  (9)  for  the  buckling  coefficient.  The  determinant  of 
equation  (51)  can  be  simplified  if  the  transverse  modulus  of  elasticity  of  the 
core  (Er)  is  assumed  to  be  infinite.  For  most  core  materials  except  possi¬ 
bly  for  low  density  foams  Ec  is  sufficiently  large  so  that  this  assumption 

yields  values  of  the  critical  pressure  that  are  only  very  slightly  too  great. 

_ 

—Underlined  numbers  in  parentheses  refer  to  the  references. 
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Before  Ec  is  allowed  to  approach  infinity,  the  first  and  fourth  columns  of  the 
-  GRz 

determinant  are  multiplied  by  — — —  and  the  third  column  is  divided  by  Ec. 


Then  when  Ec  approaches  infinity,  the  expressions  in  rows  3,  4,  5,  and  6  in 
column  3  approach  zero. 


The  expressions  in  each  row,  excepting  the  first,  are  replaced  by  new  ex¬ 
pressions,  as  indicated  by  the  following  formulas  in  which  R  represents  the 
expression  in  the  row  designated  by  its  subscript  and  in  some  column.  The 
primed  values  are  the  new  ones  to  be  substituted  for  the  old. 


i 


R 


2 


-  R2  ~  +  Rj 


^3  ~  ^3  +  ^2 

R4  =  r4  +  r2 


R5'  =  2R5  +  (n2  +  3  V  2)  R2' 

R6'  =  2R6  +  (n2  +  3~  k  2)  R2’ 


These  substitutions  cause  the  expressions  in  column  3  and  rows  2,  3,  4,  5, 
and  6  and  those  in  column  6  and  rows  3,  4,  5,  and  6  to  become  zero,  and 
the  6-by-6  determinant  is  readily  reduced  (by  minors)  to  a  4-by-4  determin¬ 
ant.  This  determinant  is  simplified  slightly  by  replacing  the  second  row  of 
expressions  by  the  second  row  minus  the  first  row. 

After  the  determinant  was  written  in  the  above  form,  a  change  in  parameters 
was  made  using  the  following  nomenclature: 

d  --  thickness  of  the  sandwich 

E  --  modulus  of  elasticity  of  facing 

^Rz  --  modulus  of  rigidity  of  core  in  the  radial  and  axial  directions 

Gr0  --  modulus  of  rigidity  of  core  in  the  radial  and  tangential  directions 

distance  between  facing  centroids 

qr  (1  ~  H2 

2Et 


h  - 

k  - 
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L  --  length  of  cylinder 

n  --  number  of  half  waves  in  the  circumference  of  the  cylinder 
q  --  the  external  critical  pressure  on  the  cylinder 
r  --  ittfcan  radius  of  the  sandwich  cylinder 


t  --  thickness  of  the  facings 


4r  (1  -  p.2)  GRe 


The  resulting  expressions  in  the  4-by-4  determinant  are: 
Row  1,  column  1 
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Row  Z,  column  1 


Row  3,  column  1 


Row  4,  column  1 


Row  1,  column  2 


Row  2,  column  2 


Row  3,  column  2 
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Row  4,  column  2 


Row  1 


Row  2 


Row  3 


Row  4 


Row  1 


i  ($  -  fr)3 

$  +  <f>  \  L  /  $2  (!+<{>) 


2(f)  +  4_V 

<£  +  <f>  3  R 


§  +  <t> 


log 


$  +  4> 


,  column  3 


±v  _i_  log  f-il 

3  R  §  +  <J>  $  +  <)> 


,  column  3 


$  +  4> 


,  column  3 


,  column  4 
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Row  2,  column  4 


Row  3,  column  4 


Row  4,  column  4 


By  setting  this  determinant  equal  to  zero  and  solving  for  k,  the  critical 
pressure  can  be  found  by: 


q  = 


2  Et 


(■  - 


k 


This  represents  a  theoretical  solution  for  the  critical  pressure  with  the 
assumption  that  the  sandwich  core  modulus  of  elasticity  is  infinite. 


The  upper  and  lower  limits  for  k  can  be  determined  as  given  in  Forest 
Products  Laboratory  Reports  1844-A  and  1844-B.  These  limits  are  for 
cylinders  with  membrane  facings  and  a  core  with  a  modulus  of  elasticity 
that  is  infinite.  The  lower  limit  is  approached  as  the  length  of  cylinders 
becomes  infinite  and  is  given  by  equation  (72)  of  Report  1844-A  as 


UL=Jl  ± 


+  ♦ 


1 

'  16 


l£ 


§  4> 


+  v 
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The  upper  limit  for  k  is  reached  as  n  becomes  infinite;  it  is 


K~  4(5  +4,)  v 

which  corresponds  to  the  critical  pressure  as  given  on  page  28  of  Report 
1844-B. 

From  this  critical  pressure  the  critical  hoop  compression  per  unit  length 
of  cylinder  is: 

N0  =  k  GR0 

which  is  the  usual  limit  imposed  on  the  edge  compression  of  sandwich  con¬ 
structions  with  membrane  facings  by  the  shear  instability  of  the  core. 


The  above  determinant  can  be  further  simplified  by  a  few  approximations 
without  any  significant  loss  of  accuracy.  Examination  of  the  parameters  § 
and  4>  shows  that  §  is  large  in  comparison  to  <f> .  Therefore,  any  terms  of 
($  +  4>)  can  be  assumed  as  equal  to  <|S  with  only  a  minor  loss  of  accuracy. 

Making  this  assumption  and  substituting  £  j  ,  the  determinant  reduces 

to  the  following  expressions: 


Row  1,  column  1 


-  —  +  2  V  1 
n^ 


Row  2,  column  1 


5  +  8Vi 


Row  3,  column  1 


n  -  1  26 

2  $ 

n  * 


n2+30  +4V0j--^  -  k 

L.  _ 
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Row  4,  column  1 


i*  [»2^eJ+  (v«Jj  -  f2  -  k 

Row  1,  column  2 

f2-')  (2k-‘)  +  i 

Row  2,  column  2 


Row  3,  column  2 


Row  1,  column  3 


2  4> 

I 


Row  2,  column  3 


Row  3,  column  3 
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Row  4,  column  3 


2  +  ^  n2  +  3  /3 
9 

Row  1,  column  4 


2  8 
n  -  Jr 


Row  2,  column  4 


-¥  * 


Row  3,  column  4 

4  h2  /3  J  -  J  ~  ^  jn2  +  3  0 

Row  4,  column  4 

*"2e[{  -!]•  ¥[»2  +  ^ 


By  setting  this  determinant  equal  to  zero,  the  approximate  solution  for  k  is: 


2  2 

To  obtain  this  value  of  k,  terms  containing  k  and  k  a  were  neglected.  It  was 
also  assumed  that  terms  (  1  ±  mo  )=  1  where  m  is  a  small  whole  number. 

Further  simplifications  can  be  made  by  neglecting  the  terms  of  +  ,  -  t-  »  -  ~  , 

8  i  j  j 

and -j  so  that  the  final  approximate  solution  for  k  is: 
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The  limits  of  this  approximate  solution  are  the  same  as  for  the  exact  solution 
if  the  terms  ±  are  assumed  equal  to  <§  . 


Description  of  Design  Curves 


Using  the  approximate  equation,  values  for  k  were  determined  for  various 

values  of  V,  n,  —  ,  and  a 2.  These  values  of  k  are  plotted  (figs.  1-5)  for 
Li  — 

values  of  —  ranging  from  1  to  100.  Curves  are  given  for  values  of  V  equal 

to  0,  0.5,  1.0,  2.0,  and  4.0.  For  each  of  these  values  of  V,  curves  are 

2  A  -4 

given  for  seven  values  of  a_  ranging  from  10-D  to  10  .  These  curves  apply 

to  sandwich  cylinders  with  equal  isotropic  facings  and  isotropic  cores  or 
orthotropic  cores  with  their  natural  axes  parallel  to  the  axial,  tangential,  and 
radial  directions  of  the  cylinders.  The  critical  pressure  is  given  by 


2  E  t 


k 


From  a  solution  of  the  exact  determinant,  it  was  found  that  the  modulus  of 
rigidity  of  the  core  associated  with  the  radial  and  axial  directions  Grz  has 


very  little  influence  on  the  critical  pressure.  It  does  not  enter  the  formulas 

Grz 

for  the  limits  of  critical  pressure.  The  ratio  R  =  — -  does  affect  the  critical 

L  GR0 

pressure  for  small  —  values  but  for  reasonable  values  of  R_the  critical  pressure 


is  affected  only  slightly.  The  ratio,  _R,  does  not  appear  in  the  approximate 
solution  for  k. 


The  accuracy  of  the  approximate  solution  was  checked  by  making  a  comparison 
of  the  approximate  and  exact  solutions  for  values  of  V  equal  to  0,  0.5,  and  1.0, 
r  t  c 

—  equal  to  50,  and  —  equal  to  1.  This  comparison  is  shown  on  the  top  curve  of 
a  a 
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figures  1,  2,  and  3  where  the  small  circles  represent  the  exact  solution.  This 

shows  that,  except  for  small  —  values,  the  error  encountered  with  the  approx- 

r 

imate  solution  is  less  than  1  percent.  For  larger  values  of  .V  and  -  the  error 

is  less  than  is  shown  here.  For  the  small  values  of  —  ,  it  is  doubtful  that 

r 

either  solution  gives  the  correct  results. 


FPL-07 


-12- 


1.0-21 


References 


(1)  Bodner,  S.  R.  ,  and  Berks,  W. 

1952.  The  Effect  of  Imperfections  on  the  Stress  in  a  Circular  Cylindrical 
Shell  under  Hydrostatic  Pressure.  PIBAL  Report  No.  210, 
illus.  Polytechnic  Institute  of  Brooklyn. 

(2)  Donnell,  L.  H. 

1956.  Effect  of  Imperfections  on  Buckling  of  Thin  Cylinders  under  Exter¬ 
nal  Pressure.  Jour.  Appl.  Mech.  ,  23:  (No.  4). 

(3)  Galletly,  G.  D. ,  and  Bart,  R. 

1956.  Effects  of  Boundary  Conditions  and  Initial  Out-of-RoundneBs  on 

the  Strength  of  Thin-Walled  Cylinders  Subject  to  External  Hydro¬ 
static  Pressure.  Jour.  Appl.  Mech.,  23:  (No.  3). 

(4)  Kempner,  Joseph,  Pandolai,  K.  A.  V.  ,  Patel,  S.  A.,  and  Crouzet-Pascal,  J. 

1957.  Post-Buckling  Behavior  of  Circular  Cylindrical  Shells  Under  Hydro¬ 

static  Pressure.  Jour.  Aero.  Sci.  ,  24:  (No.  4). 

(5)  Kirstein,  A.  F. ,  and  Wenk,  E. 

1956.  Observations  of  Snap-Through  Action  in  Thin  Cylindrical  Shells 
under  External  Pressure.  David  Taylor  Model  Basin  Report 
1062.  U.  S.  Navy  Dept, 

(6)  Langhaar,  H.  L.  ,  and  Boresi,  A.  P. 

1955.  Snap-Through  and  Post-Buckling  Behavior  of  Cylindrical  Shells 

under  the  Action  of  External  Pressure.  T.  &  A.M.  Report  No. 

80.  Department  of  Theoretical  and  Applied  Mechanics,  Univer¬ 
sity  of  Illinois. 

(7)  Mises,  R.  von. 

1929.  The  Critical  External  Pressure  of  Cylindrical  Tubes  under  Uniform 
Radial  and  Axial  Load.  Stodola-Festschrift,  Zurich. 

(8)  Nash,  William  A. 

1955.  Effect  of  Large  Deflections  and  Initial  Imperfections  on  the  Buck¬ 
ling  of  Cylindrical  Shells  Subject  to  Hydrostatic  Pressure.  Jour. 
Aero.  Sci.,  22:  (No.  4). 

(9)  Raville,  Milton  E. 

1954.  Analysis  of  Long  Cylinders  of  Sandwich  Construction  under  Uniform 
External  Lateral  Pressure.  U.S.  Forest  Products  Laboratory 
Reports  No.  1844,  23  pp. ,  illus.;  and  supplementary  Reports  No. 
1844-A  (1955)  and  1844-B  (1955). 


FPL-07 


-13- 


(10)  Timoshenko,  S, 

1936.  Theory  of  Elastic  Stability.  McGraw-Hill,  New  York. 

(11)  U.  S.  Forest  Products  Laboratory 

1955.  Sandwich  Construction  for  Aircraft.  Part  H,  Second  Edition  1955. 

ANC-23  Bulletin,  115  pp. ,  illus.  Superintendent  of  Documents, 
U.  S.  Government  Printing  Office,  Washington  25,  D.  C. 

(12)  Volmir,  A.  S. 

1957.  On  the  Influence  of  Initial  Imperfections  on  the  Stability  of 

Cylindrical  Shells  under  External  Pressure.  Dokladi  Akad. 
NaukSSSR  (N.  S. )  113:  (2)291-293  (translated  by  M.  D. 
Friedman,  Inc.,  Needham,  Mass.,  V-124,  5  pp.) 

(13)  Wenk,  E. ,  Slankard,  R.  C. ,  and  Nash,  W.  A. 

1954.  Experimental  Analysis  of  the  Buckling  of  Cylindrical  Shells  Sub¬ 
jected  to  External  Hydrostatic  Pressure.  Proc.  Soc.  Exp. 
Stress  Anal.,  XII  (No.  l). 

(14)  Windenburg,  D.  F. ,  and  Trilling,  C. 

1934.  Collapse  by  Instability  of  Thin  Cylindrical  Shells  under  External 
Pressure.  Trans.  American  Society  of  Mechanical  Engineers, 
56  (No.  11). 


FPL-07 


-14- 


FPL -07 


Figure  2. — Values  of  k  for  V  =  0  5. 
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Figure  3. --Values  of  k  for  V  =  1.0. 


Figure  4, --Values  of  k  for  V  =  2.0. 
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Figure  5. --Values  of  k  for  V  =  4.0. 
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